Estimation of a classical parameter with gaussian probes: magnetometry with 

collective atomic spins 
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We present a theory for the estimation of a classical magnetic field by an atomic sample with a 
gaussian distribution of collective spin components. By incorporating the magnetic field and the 
probing laser field as quantum variables with gaussian distributions on equal footing with the atoms, 
we obtain a very versatile description which is readily adapted to include probing with squeezed 
light, dissipation and loss and additional measurement capabilities on the atomic system. 
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External classical perturbations of a quantum system 
cause changes in the state of the system, and a measure- 
ment of a suitable observable provides an estimate of the 
strength of the perturbation. Atoms are excellent probes 
for the estimation of, e.g., classical electric and magnetic 
fields as well as for rotations and accelerations of inertial 
frames. The formal description of such ultra-sensitive 
measurements is quite complicated and has only been 
formulated recently. The main difficulty arises from the 
fact that the quantum state of the atoms is changed due 
to both the interaction with the classical perturbation 
and the measurement process itself which yields a time 
series of stochastic outcomes. Quantum trajectory the- 
ory Q makes it possible to simulate this stochastic pro- 
cess, and descriptions are available which combine the 
quantum dynamics and the parameter estimation con- 
ditioned on the detection record j3, Recently, the 
classical theory of Kalman filters was combined with the 
quantum trajectory theory and under the assump- 

tion that the quantum state of the atomic system could 
be treated as a gaussian state of oscillator-like degrees 
of freedom, and the initial uncertainty about an applied 
magnetic field could also be described by a gaussian dis- 
tribution function, analytical expressions for the preci- 
sion of the estimate of the field were derived. The analy- 
sis showed that the probing of the atomic system squeezes 
the atomic observable and results in a measurement un- 
certainty that decreases with time t and atomic number 
N^t as l/(-/Vati'^^^) and not as l/V^at^, as one might have 
expected from standard counting statistics arguments. 

Here, we present an alternative quantum theory for the 
estimation of a i?-field by an atomic probe. The idea is 
to treat both the laser field used to probe the atoms, the 
atoms themselves, and the classical i?-field as one large 
quantum system. Quantum mechanical state reduction 
associated with measurements then provides directly the 
estimate for the expectation value and uncertainty for the 
quantity of interest. Our theory arrives easily at final 
estimation results, and it readily generaHzes to include 
decay and losses. 



We will assume that a gaussian state, fully character- 
ized by expectation values and covariances, describes the 
laser field, the atoms and the i?-field, and we will use that 
the gaussian character of the state is preserved during the 
evolution due to the interactions and measurements in- 
volved. We benefit from the considerable attention given 
to the transformation of gaussian states under interac- 
tions and measurements because this class of states per- 
mits a detailed characterization of entanglement issues 
(see, e.g., and references therein). 

We consider a collection of atoms with a spin-1/2 
ground state, polarized along the a;-axis. The i?-field 
is assumed to point along the y-axis, and it hence causes 
a Larmor rotation of the atomic spins towards the z-axis. 
A linearly polarized optical probe is transmitted through 
the gas. The linear probe is decomposed into two cir- 
cular components, and different couplings to an excited 
state introduce a phase difference of the two field com- 
ponents and cause a Faraday rotation of the polariza- 
tion proportional to the population difference between 
the atomic ground states. It is the recording of 
this rotation that enables us to determine the B-field. 
The atoms are effectively described by a collective spin 
operator J = f X^i'''^'^ ^nd the polarization compo- 
nents of the field are described by a Stokes vector S. 
With the initially spin-polarized sample, and the inci- 
dent field in a linearly polarized state, we may treat 
Jx and Sx as classical variables related to the number 
of atoms iVat and photons iVph via {Jx) = -^^^ and 



{Sx) — ^^^'^ . When the field is not too close to reso- 
nance, we may eliminate the excited states, and the ef- 
fective Hamiltonian of the atom-light interaction can be 
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written as H oc 2^ J^Sz, with A the detuning from res- 
onance. The coupling strength between a single atom 
and the radiation field (quantized within a segment of 

length L = CT and area A) is g = /^^^ d/ H with d 
the atomic dipole moment and hio the photon energy. It 
is convenient to introduce effective dimensionless posi- 
tion and momentum operators for the non-classical com- 
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ponents of the spin and Stokes vector, Xat 

s. 



Pat 



Jz 



X 



tators [xi,a;j] = \pi,Pj\ = 0, [a;i,pj] = i5ij. The perfectly 
polarized atomic state and the laser field polarized along 
the a;-direction correspond to the ground state, i.e., a 
gaussian minimum uncertainty state of the harmonic os- 
cillator associated with these variable. 

We assume that the probing of the atoms takes place 
with a continuous wave field. Such a field can be treated 
as a succession of beam segments of duration r and with 
a given mean number of photons iVph — 2{Sx) /h—^Tin 
each segment, with <& the photon fiux. The continuous 
measurement of the field is then broken down into indi- 
vidual measurements on each segment. The continuous 
limit is achieved when r ^ and Nph in each segment 
gets correspondingly small. In the limit of small r, the 
integral over r is equivalent to the application of a coarse 
grained Hamiltonian given by Ht = ftwrPatPph with di- 



with commu- 



mensionless = ^\/^^r = i<i>^3/2_ 
Due to the r-dependence of g, Kr is proportional to y/r. 
When we incorporate the S-field coupling to the atoms, 
PBJy/h, with l3 the atomic magnetic moment, the total 
effective Hamiltonian is given by 



Ht = h {KrPatP ph + Mr-BXat) , 



(1) 



^ith = i/3y^T. 

We treat the classical B-field variable on equal footing 
with the quantum variables. The Heisenberg equations 
of motion for the column vector of the five variables y = 
(S,Xat,Pat,a;ph,Pph)^ yield y{t + t) ^ Sry{t) with the 
transformation matrix 

/I \ 

1 

Sr= -Air 1 . (2) 

Kr 1 

\ 1 / 

The covariance matrix, defined as in 0, Q, 7y — 



2Re{{yi — {yi)){yj — (yj))) then transforms as 

Jit + t) = Srj{t)Sj , 



(3) 



due to the atom-light and the atom-field interaction. In 
the gaussian approximation, the system is fully charac- 
terized by the vector of expectation values (y) and the 
covariance matrix 7. We probe the system by measuring 
the Faraday rotation of the probe field, i.e., by measur- 
ing the field observable a;ph. Since the photon field is an 
integral part of the quantum system, this measurement 
will change the state of the whole system, and in particu- 
lar the covariance matrix of the residual system of atoms 
and B-field. We denote the covariance matrix by 
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(4) 



where the 3x3 sub-matrix is the covariance matrix 
for the variables yi = (B,Xe,t,Pat)'^ , is the 2x2 
covariance matrix for 1/2 = {xph,Pph)^ , and is the 
3x2 correlation matrix for yi and t/l". An instantaneous 



measurement of 



^ph 



then transforms A^ as 



A' = A^ - C^(7rB-y7r) 



(5) 



where tt = diag(l,0), and where the inverse denotes the 
Moore-Penrose pseudoinverse, as (ttB^^tt) is not invert- 
ible. Equation Q is equivalent to the result for classical 
gaussian random variables derived, e.g., in After the 
measurement, the field part has disappeared, and a new 
beam segment is incident on the atoms. This part of the 
beam is not yet correlated with the atoms, and it is in 
the oscillator ground state, hence the covariance matrix 
7 is updated with A'^, C^^ a 3 x 2 matrix of zeros, and 
B; = diag(l,l). 

Unlike the covariance matrix update, which is indepen- 
dent of the value actually measured in the optical detec- 
tion, the vector (y) of expectation values will change in 
a stochastic manner depending on the outcome of these 
measurements. The outcome of the measurement on Xph 
after the interaction with the atoms is random, and the 
actual measurement changes the expectation value of all 
other observables due to the correlations represented by 
the covariance matrix. Let x denote the difference be- 
tween the measurement outcome and the expectation 
value of Xph, i.e., a gaussian random variable with mean 
value zero and variance 1/2. The change of (yi) due to 
the measurement is now given by: 



{yi) 



{Vi) = 



= (yi) + C,(7rB7r)-i(x,0)^ 



(6) 



where we use that the measurement on Xph only, 
leads to the particularly simple form (7ri?7r)~^ = 
diag(i3(l, 1)~^,0), and hence the actual value of the sec- 
ond entrance in the vector (X7O) is unimportant. 

The gaussian state of the system is propagated in 
time by repeated use of lO and the measurement up- 
date formulae (jHJ-®- This evolution is readily imple- 
mented numerically, and the expectation value and our 
uncertainty about the value of the i?-field are given by 
the first entrance in the vector of expectation values 
(yi) — (B) and the (1,1) entrance in the covariance ma- 
trix Aj{l, 1) = 2(AB)2. 

The above discussion specifies how the parameter 
estimation can be performed. In the problem at 
hand, the variable x^t does not couple to B and 
Pat, and we are left with a closed 2x2 system for 
the reduced covariance matrix of B and p^t- V = 
[2(AB)2,2(ASpat)';2(Apat5)',2(Apat)']- In the limit 
of infinitesimally small steps the update formulae ©-© 
translate into a differential equation on the matrix Ri- 
catti form 

V(t) = -DV(t) - V(t)D^ - V(<)EV(i), (7) 
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Figure 1: Uncertainty of the B-field in pT (IpT = 10"^^ T) 
as function of time. We use a 2 mm^ interaction area, 2 x 10^^ 
atoms, 5 x 10^^ photons s~^, ABq = 1 pT, 1 GHz detuning, 
and 852 nm light, appropriate for the ^^^Cs(6S'i/2 (-F = 4) - 
GPi/2iF = 5)) transition with decay rate 3.1x10^ and cor- 
responding atomic dipole moment d = 2.61 x 10~^^ Cm. The 
effective couplings are « = 1.83x10'' s"^ and fi = 8.79x 10* (s 
pT)~^. Factors of order unity related to the coupling matrix 
elements among different states of the actual Zeeman sub- 
structure are omitted. The lower curves are without inclusion 
of atomic decay, and the upper curves include atomic sponta- 
neous emission with a rate rj — 1.7577 s~^. The solid (dashed) 
curves are for coherent (squeezed, r — S) optical probe fields 
(see text). 



with D = [0, 0;^,0], E = diag(0,K^), k = k^/v^, and 
/i = Mt/t- We solve Q by expressing it in terms of 
two coupled linear matrix equations W = — DW, U = 
EW + D'^U, V = WU^i [H, and find the analytical 
solution for the variance of the magnetic field 



AB{tf 



{1 + kH)AB^ 



(8) 

with ABq the initial variance. In the limit of K^t ^ 1, 
we have AB{ty ~ 6/(^2^2^3^ explicitly giving the l/7Vat 
and l/t~^ scaling also found in 

The lower solid curve in Fig. 1 shows the uncertainty 
of the i?-field as a function of time. It is worth pointing 
out that compared with 0|, not only the spirit in which 
we deal with J5 as a quantum variable but also the formal 
derivation is different. In the Kalman filter equation 
deals with the covariance matrix for the joint estimator 
of the classical B-field and the mean value of the atomic 
spin component along the z-axis. The latter variance is 
initially zero, because we assume that the mean value is 
initially known to be zero. Our covariance matrix deals 
with two quantum observables, and neither have a van- 
ishing variance in the initial state. 



We may now go back to (EJ and derive the stochastic 
differential equation 

d{Bit)) = V2K{ABp^tfdW{t) (9) 

for the expectation value of the i?-field. Here dW{t) = 
xV^di is a Wiener increment with gaussian white-noise 
statistics {dW{t)) = 0, (dWitY) = dt. {ABp^^? =^ 
3/{K^^t^) in the long time limit as determined by the 
Ricatti equation 0) ^^'^ it follows that the locking of 
the value of (B), conditioned on the measurements, takes 
place predominantly in the early stages of the detection 
process. This is in agreement, of course, with the rapid 
reduction of the uncertainty as a function of time. 

Together with the phase shift, there is a small proba- 
bility that the atoms decay by spontaneous emission from 
the upper probe level to one of the two ground states. 
This occurs with a rate 77 = <I>-J ^ ryJ+A^ ) ' '^here F is 

the atomic decay width and a = )? /{2tt) is the resonant 
photon absorption cross-section. The consequence of the 
decay is a loss of spin polarization. If every atom has 
a probability 77,- = rjT to decay in time r with equal 
probability into the two ground states, the collective 
mean spin vector is reduced by the corresponding fac- 
tor (J) («/)(! — r]r). When the classical x-component 
is reduced this leads to a reduction with time of the cou- 
pling strengths Kr ^ K,r\f^ 
which was also discussed in 
pectation values evolves as (y(t 
= diag(l, VI - ^r, VI - '^T, 1, 1). 
The fraction r\r of atoms that have decayed repre- 
sents a loss of collective squeezing because its correlation 
with the other atoms is lost, whereas it still provides a 
contribution Y? jA, per atom to the collective spin vari- 
ance. The mean value of J2 can be expressed in terms 
of the mean values of the iVat(A^at — 1) atomic correla- 
tions a\a\, and counting terms, we find that (J^) 
{\-^rflyJl) + (ft2Arat/4)(l- (1 -?7^)2). Translating this 
and similar expressions for J2 a^d Jx-Jy into the appro- 
priate formulae for the effective position and momentum 
observables, (jnj generalizes to 



and /ir 1-^ /ir-v/l — ??r , 
, and the vector of ex- 
r)) =L,S,(y(i)) with 



7(t + r) =L,S,7(i)S^L, 



(10) 



for r\r with M.^ — diag(0, r\r,"r\T-, 0, 0). The prefactor 
f^N&t/ {Jx{t)) initially attains the value 2, and increases 
by the factor (1 — 1]^-)^^ in each time step r. The effects 
of measurements on the covariance matrix and the expec- 
tation value vector are obtained as in the case without 
noise, and for ry^ = we regain the noise-less case. 

The upper solid curve in Fig. 1 shows the results of 
the measurement when noise is taken into account. The 
covariance matrix makes the atomic probe broader, and 
simultaneously, the effective coupling of the atoms to the 
light field and to the i?-field is reduced, so that the knowl- 
edge acquired in the initial detection stages is preserved 
but the uncertainty AB does not decrease indefinitely. 
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The value of B is estimated by the polarization ro- 
tation of the optical field, and it is natural to enquire 
whether the use of polarization squeezed Hght with a 
smaller variance of Xpi^ may be utilized to improve the 
estimate. To analyze this proposal, we go back to our 
update formulae and represent each new segment of the 
incident field with gaussian variances B^^ — diag(l/r, r), 
and leave all other operations unchanged. The result is 
a reduction of the variance of our estimate, shown as the 
dashed curves in Fig. 1. The upper dashed curve is for 
the case when noise is included. The Ricatti equation 
can be solved in the noise-less case, and the only change 
of the result in l|Hl is that all occurrences of are re- 
placed by tk'^. In the long time Hmit, the estimate is 
improved by the factor 1/r. Since the optical field is 
not squeezed if the time segments r are shorter than the 
squeezing bandwidth fl, we rely on a separation of time 
scales fl^^ ^ T ^ M^^7 1^ '^ for the above update formu- 
lae to be vaHd, and for the Ricatti equation to provide 
a precise analytical solution. For the parameters used in 
Fig. 1, the squeezing bandwidth should be larger than 
10 MHz. Effects of finite squeezing bandwidth will be 
analyzed elsewhere. 

We can improve our estimate by noting that the covari- 
ance matrix describes correlations between the atomic 
observables and the i3-field, and the uncertainty in the 
measurement is linked with the uncertainty of the atomic 
observable x^t- After the optical probing it is in prin- 
ciple possible to perform a destructive (Stern-Gerlach) 
measurement of this atomic variable. This can of course 
only be done once. The formal treatment of measure- 
ments in also applies when the atomic component is 
being measured, and we can readily determine the new 
variance on the i?-field estimate. From the Ricatti equa- 
tions we know the covariance matrix A-y analytically, and 
assuming an atomic measurement at time t, we obtain 

^BsG{t)^ = 1+2^2 (ABo)^^^+§Iv^(Ai3o)^«^ - This variance 
is smaller than (A_B(t))^ from and in the long-time 
limit the variance is reduced by a factor of 4. 

In summary, we have described a theory for the es- 
timation of a classical i?-field by an atomic ensemble 
with a gaussian distribution of collective spin compo- 
nents. Our theory makes use of results obtained in the 
study of classification and characterization of entangle- 
ment in continuous- variable systems In general, the 
gaussian ansatz holds for Hamiltonians which are at most 
second order polynomials in the canonical variables, and 
the gaussian character of a system is maintained under 
physical operations which are implemented using linear 
optical elements and homodyne measurements It is 
clearly convenient to have a unified formalism that deals 
with both the probing field, the atomic probe, and the un- 
known B-field, and which bypasses the need for separate 
probabilistic arguments to yield the final estimator. The 



treatment of the unknown -B-field as a quantum variable 
is not incompatible with our assumption that it is a classi- 
cal parameter. We may imagine a canonically conjugate 
variable to B having an uncertainty much larger than 
required by Heisenberg's uncertainty relation and/or ad- 
ditional physical systems, entangled with the B-variable, 
in which cases the i?-distribution is indeed incoherent 
and "classical". Also, one may argue that all classical 
variables are actually quantum mechanical variables for 
which a classical description suffices, and hence our the- 
ory provides the correct estimator: quantum mechanics 
dictates that the quantum state provides all the available 
knowledge about a system, and any estimator providing 
a tighter bound hence represents additional knowledge 
equivalent to a hidden variable, and this is excluded by 
quantum theory. It is of course crucial that our measure- 
ment scheme corresponds to a quantum non-demolition 
(QND) measurement, i.e., we assume that there is not 
a free evolution of the i?-field induced by its conjugate 
variable which may thus remain unspecified. It is also 
this QND property of the measurement scheme that im- 
plies the monotonic reduction of AB which is consistent 
with the classical parameter estimation (we can not un- 
learn what we have already learnt about B), unHke, e.g., 
the uncertainty of the atomic x^t variable which must in- 
crease when Apat is reduced and when the atoms undergo 
spontaneous decay. 

We expect extensions of the present theory to be appli- 
cable to the description of a variety of experiments aim- 
ing at ultra-high precision, including, e.g., atomic clocks, 
studies of parity violation, and the detection of gravita- 
tional waves. 
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